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ABSTRACT 

With the excellent angular resolution of the Chandra X-ray Observatory, it is possible to geomet- 
rically determine the distance to variable Galactic sources, based on the phenomenon that scattered 
radiation appearing in the X-ray halo has to travel along a slightly longer path than the direct, un- 
scattered radiation. By measuring the delayed variability, constraints on the source distance can be 
obtained if the halo brightness is large enough to dominate the point spread function (PSF) and to 
provide sufficient statistics. The distance to Cyg X-3, which has a quasi-sinusoidal light curve, has 
been obtained with this approach by Predehl et al. Here we examine the feasibility of using the 
delayed signature of type I X-ray bursts as distance indicators. We use simulations of delayed X-ray 
burst light curves in the halo to find that the optimal annular region and energy band for a distance 
measurement with a grating observation is roughly 10-50" and 1-5 keV respectively, assuming Chan- 
dra's effective area and PSF, uniformly distributed dust, the input spectrum and optical depth to 
GX 13+1, and the Weingartner & Drainc interstellar grain model. We find that the statistics are 
dominated by Poisson noise rather than systematic uncertainties, e.g., the PSF contribution to the 
halo. Using Chandra, a distance measurement to such a source at 4 (8) kpc could be made to about 
23% (30%) accuracy with a single burst with 68% confidence. By stacking many bursts, a reasonable 
estimate of systematic errors limit the distance measurement to about 10% accuracy. 
Subject headings: X-rays: ISM — X-rays: bursts — dust, extinction — stars: individual (GX 13+1) 



1. INTRODUCTION 

Type I X-ray bursts potentially provide a delayed sig- 
nature in the halo that can be used to measure the dis- 
tance to variable X-ray sources with a method proposed 
by Trumper & Schonfelder (1973; hereafter TS73), prior 
to the detection of the first X-ray halo around GX 339- 
4 by Einstein (Rolf 1983). The TS73 method is based 
on the phenomenon that grain-scattered photons from 
an X-ray source travel along slightly longer paths to the 
telescope than the direct, unscattercd photons. Tempo- 
ral variations in the intensity are therefore delayed and 
smeared when they appear in the X-ray halo. Using the 
single scattering approximation, it is possible to predict 
the delayed light curve as a function of halo angle, energy, 
and distance to an accuracy that is limited by knowledge 
of the distribution of dust along the line of sight and by 
the accuracy of the interstellar grain model. Hu et al. 
(2004) developed a similar method to measure distances 
with X-ray halo variability using the frequency domain. 

The feasibility of the TS73 distance measurement 
method has been demonstrated by Predehl et al. (2000), 
who derived a distance of 9^2 kpc to Cyg X-3 using the 
source's quasi-sinusoidal variability. In this paper, we ex- 
amine the applicability of the TS73 method to the abrupt 
increase in flux provided by type I X-ray bursts, which 
result from unstable thermonuclear ignition of accreted 
material on the surface of the neutron star. Sharp and 
temporary intensity changes, such as that produced by 
type I bursts or gamma-ray bursts (GRBs), produce a 
qualitatively different response in the X-ray halo than 
the response from smooth intensity variations like those 
of Cyg X-3. For example, consider a smooth change in 
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point source flux over an hour or so. Eventually the 
halo brightness will reach a new base level correspond- 
ing to the new point source intensity, and information 
on the source distance and dust distribution is contained 
in the rate of response and the resulting phase shift of 
the light curve. With X-ray bursts a new base level 
is not obtained; rather, the burst appears as a slight 
and temporary increase in halo flux after a time delay 
that depends on the scattering geometry. The intrinsi- 
cally small signal-to-noise (S/N) ratio provided by bursts 
makes their potential application as distance indicators 
particularly difficult, requiring the physical parameters 
governing the delayed light curve in the X-ray halo be 
measured to a high degree of accuracy, and requiring 
careful treatment of the point spread function (PSF). 

Gamma-ray bursts produce a similar halo light curve, 
although the problem is much simpler: With spectro- 
scopic determination of the distances to host galaxies of 
GRBs, the degeneracy between the distribution of dust 
and the distance to the source is removed. This has been 
exploited by both the Swift and XMM-Newton observa- 
tories to determine the distances to Galactic interstellar 
dust clouds along various sightlines (Vaughan et al. 2004; 
Tiengo & Mereghetti 2006; Vaughan et al. 2006). 

The physical parameters governing the delayed light 
curve in the halo are the peak and persistent flux, the 
spectral shape, the scattering optical depth, the scale 
time of the approximately exponential burst light curve, 
and the distribution and composition of interstellar dust 
along the line of sight. Two of the most widely used in- 
terstellar grain models are the Mathis, Rumpl, & Nord- 
sieck (1977; MRN) model and the Weingartner & Draine 
(2001; WD01) model. The MRN model is composed 
of silicate and graphite grains with a size distribution 
of n(a) oc a~ 3 5 , which reproduces the observed extinc- 
tion of starlight. The WD01 grain model additionally 
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accounts for the observed infrared and microwave emis- 
sion from the diffuse ISM by including sufficient small 
carbonaceous grains. In this paper we use the WD01 
interstellar grain model and we assume that the dust is 
distributed uniformly along the line of sight. 

The structure of this paper is as follows: We begin by 
describing the Chandra observation that was used in this 
analysis and we briefly discuss pile-up (§2). In § 3 we 
present two analytic equations from Draine (2003) de- 
scribing the scattering of X-rays from interstellar grains 
and the resulting halos, from which we derive an equa- 
tion describing the predicted counts spectrum for an ar- 
bitrary range of angles and energies. We also discuss our 
treatment of the PSF. The derived equation is used to 
estimate the optical depth to GX 13+1 assuming uni- 
formly distributed dust along the line of sight by fitting 
it to the observed surface brightness distribution (SBD). 
In § 4 we demonstrate our method for calculating the ex- 
pected delayed halo light curve of singly-scattered pho- 
tons from a burst. In § 5 we show that, given a particular 
set of assumptions, the 1-5 keV energy range and the 10- 
50" annular region yield the highest time-averaged S/N 
for viewing bursts in the halo. With this binning, simu- 
lations are carried out to test the feasibility of applying 
the TS73 method to type I X-ray bursts. We finish in 
§ 6 by summarizing our results. 

2. OBSERVATIONS 

In this paper, we use a 29 ks Chandra observation of 
GX 13+1 (Obs. ID 2708). The observation uses the 
Advanced CCD Imaging Spectrometer (ACIS), with the 
source focused on the back-illuminated ACIS-S3 chip, 
and with the High Energy Transmission Grating (HETG) 
inserted in the optical path. GX 13+1 has a substantial 
halo (Smith et al. 2002), yet it is not a consistent type 
I burster (two bursts were observed in Sept. 1989; Mat- 
suba et al. 1995). Nevertheless, this observation allows 
us to test the potential for using bursts as distance indi- 
cators. We accomplish this by measuring the GX 13+1 
SBD, from which we predict the size of a delayed burst 
signal in the halo as a function of delay time. Figure 

shows the observed surface brightness distribution of 
GX 13+1 (diamonds with error bars) and the best-fit 
halo (see § 3). Statistical errors for the SBD are slightly 
less than 1% for 9^ > 10". Included in the figure is the 
contribution from the X-ray background (XRB) and the 
PSF, the treatment of which is in § 3. Note the clear 
turnover in the observed SBD within the innermost 2- 
3", resulting from the rejection of events due to pile-up. 
Pile-up is defined to be when more than one photon ar- 
rives in the same 3x3 pixel island in a single frame time. 
This can result in spectral hardening, "grade migration," 
or rejected events in the case of pulse saturation. Accord- 
ing to the Chandra Proposers' Observatory Guide v8.0 
(CPOG), a pile-up fraction 2 of 10% will very likely im- 
pact the observation. Smith et al. (2002) found the pile- 
up fraction reaches 10% when the count rate per pixel 
is 0.012 counts/frametime/pixel, which corresponds to 
halo angles less than about 3" in this observation. While 
pile-up in the GX 13+1 observation is most prevalent for 
halo angles < 5", it likely affects the data to larger angles. 
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Fig. 1. — Surface brightness distribution of GX 13+1 (di- 
amonds with error bars). The halo profile assumes uniformly 
distributed dust. The best-fit column density for 10-100" is 
JVh = (2.75±0.06)xl0 22 cm -2 , corresponding to T sca (l keV) w 1.4 
using the WD01 grain model. For halo angles less than 10" the 
PSF is modeled using ChaRT, and for halo angles greater than 
10" the PSF is modeled using a power law with an exponential 
term, using coefficients from Gaetz (2004). The pile-up fraction 
reaches 10% at about 0.012 counts/frametime/pix (Smith et al. 
2002), corresponding to 0^ < 3". The lower panel shows the ratio 
of the halo flux to the sum of the PSF and background fluxes (solid 
curve), and the fractional difference in the halo contribution to the 
SBD assuming half of the dust is uniformly distributed, and half is 
divided into two clumps at x = 0.33 and x = 0.66 (dashed curve). 



Given the relatively large calibration uncertainties of the 
PSF (CPOG) and its increasing contribution nearer to 
the point source, we conservatively ignored halo angles 
less than 10" in our final analysis. 

Spectral response matrices and auxiliary response files 
for GX 13+1 were created with the standard CIAO 
scripts. 3 The GX 13+1 source spectrum was fitted using 
XSPEC vll.3.2 with an absorbed blackbody plus a power 
law model, which allowed us to infer its hydrogen column 
density, measure the source flux, and predict X-ray halo 
spectra. The best-fit spectral parameters for this obser- 
vation are kT hh = 0.81 ± 0.02 keV, T = 0.81±gJ|, and 
A H = (2.46 ± 0.04) x 10 22 cm~ 2 . 

3. GRAIN-SCATTERED X-RAY HALOS AND THE PSF 

The theory of X-ray scattering and the resultant X- 
ray halos has been presented in detail in many previous 
publications (e.g., Overbeck 1965; Mauche & Gorenstein 
1986; Klose 1991; Mathis & Lee 1991; Smith & Dwek 
1998; Draine 2003). In this analysis, we use two ana- 
lytic equations of Draine (2003, eqs. [1] and [2] below), 
we assume uniformly distributed dust, and we apply the 
WD01 interstellar grain model 4 with Ry — 3.1. Using 
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the WD01 grain model, Draine (2003) obtained an ap- 
proximation to the differential scattering cross section 
that applies to photon energies greater than 0.5 keV: 



da(<t>,E) 
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where SC a,5o ~ 360"/i?(keV) is the median scattering 
angle, which encloses half of the halo flux. This type 
of dependence on scattering angle <j) implies that higher 
energy photons lead to halos that are more centrally 
peaked, and also implies that a greater fraction of the X- 
ray halo at higher photon energies results from smaller 
angle scatterings occurring close to the observer. This 
has observable implications in the delayed light curve 
following a type I burst (§ 4). The fraction of the total 
single-scattered halo flux interior to halo angle #h is char- 
acterized by the normalized (dimensionless) dust density 
p(x): 
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where x is the fractional distance to the source. Using 
these relations, we define the fractional halo density (in 
units of arcsec -2 ) as 
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With an appropriate normalization, the fractional halo 
density for a given energy will produce approximately the 
same halo profile as J sca (fl), the conventional formula for 
halo surface brightness distributions (e.g., Mathis & Lee 
1991). With its explicit energy and halo angle depen- 
dence, however, this new quantity is particularly useful 
for calculating the delayed burst light curve in the halo. 
The spectrum of the X-ray halo integrated over all halo 
angles is a function of the source's spectrum prior to ab- 
sorption and scattering along the line of sight F*(E), the 
measured point source spectrum Fx {E) , and the optical 
depths through: 
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where T abs ( sca ) is the optical depth to absorption (scat- 
tering) and F±e~ T * ba — Fx + -Fhaio- As discussed below, 
the distribution function for delayed photons scattering 
along the line of sight has a simple analytic form only for 
singly-scattered photons, yet the probability of a photon 
emerging in the halo after undergoing exactly N scatter- 
ings is r^ a /N\ relative to the observed source strength, 
so multiple scattering is likely important for sources with 
appreciable halos (Mathis & Lee 1991). Nevertheless, 
singly-scattered photons generally appear earlier in the 
halo (after a shorter time delay) than photons that have 
been scattered multiple times (see Fig. 9, Alcock & 
Hatchett 1978), so the effects of multiple scattering can 
be minimized by choosing small delay times. The halo 
spectrum for singly-scattered photons is obtained by re- 
placing the term [e Taca — 1] on the right side of eq. (0J 
with r sca . Figure |2Ja shows the measured point source 
spectrum, the inferred halo spectrum, and the single- 
scattered halo spectrum for GX 13+1. For an arbitrary 
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Fig. 2. — (a) Best-fit spectral model for GX 13+1 (bold solid 
curve), the halo spectrum integrated over all halo angles (solid 
curve), and the single-scattered halo spectrum (dotted curve), using 
the WD01 interstellar grain model, (b)-(c) The predicted counts 
spectra, obtained by multiplying the halo spectra by the Chandra 
zeroth-order effective areas for grating and non-grating observa- 
tions. The decrease in effective area with off-axis angle has been 
ignored. For 0^ < 100" and E < 5 keV, the decrease in area is 
only ~1%. 



energy range, the SBD as observed by Chandra (minus 
the PSF contribution) can be obtained from 
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By removing the fractional halo density from eq. (|SJ), 
equivalent to integrating over all halo angles, the inte- 
grated counts spectrum is obtained (Fig. 26 and Fig. 
2c). 

To accurately model the X-ray halo, it is extremely 
important to model the PSF so that the relative contri- 
butions of the dust-scattered and mirror-scattered halos 
can be separated. The Chandra Ray Tracer 5 (ChaRT) 
can be used to obtain a ray-traced PSF for any point on 
the detector and can be tailored specifically to the source 
spectrum; but while the simulated PSF represents the 
core of the source region well, it systematically under- 
estimates the PSF at larger halo angles. We therefore 
chose to use the ray-traced PSF for halo angles less than 
10", normalizing the ray-traced PSF to the GX 13+1 
source flux in 100 eV bands. Following Gaetz (2004), the 
Chandra PSF flux for halo angles greater than 10" can be 
approximated by a power law with an exponential cutoff. 
Over a given energy range, the PSF can be approximated 

5 see |http: / /cxc. harvard.edu/chart / threads/index.html| 
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by 
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where /eef.io^) is the encircled energy fraction at 
10" and is given in the CPOG for various energies, and 
the denominator is defined as 
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where the units of angle are arcsec. The denominator 
of eq. is simply a normalization; -gfgP is the inflec- 
tion point of the power law times exponential function. 
The coefficients j(E) and c(E) were measured in 1 keV 
bands by Gaetz (2004) using observations of Her X-l. 
Integration was performed by interpolating the values of 
/eef,io(-E)i j{E), and c(E) to form a nearly continuous 
grid. 

The scattering optical depth to GX 13+1 was found by 
matching the observed SBD to that produced by the sum 
of eqs. (J5J and © and the measured XRB (see Fig. 0). 
Note, however, that the use of eq. (0 assumes photons 
scattered two or more times produce the same halo SBD 
as singly-scattered photons, although multiply-scattered 
photons tend to produce somewhat broader halo profiles 
(Mathis & Lee 1991). The best-fit hydrogen column den- 
sity leading to the observed halo is (2.75 ± 0.06) x 10 22 
cm~ 2 (assuming a 20% [10%] uncertainty in the PSF 
[XRB] contribution), corresponding to T sca (l kcV) « 1.4 
using the WD01 interstellar grain model, while the best- 
fit column density using the a blackbody plus power law 
model is 7V H = (2.46 ± 0.04) x 10 22 cm- 2 (§ 2). The 
8-14% difference could be attributed to an incorrect as- 
sumption of uniformly distributed dust, to an inaccurate 
grain model, to the assumption that the SBDs for singly- 
and multiply-scattered photons are the same, or perhaps 
the gas-to-dust ratio along the sightline differs from the 
Galactic mean. 

Finally, the measured halo SBD can also be used to 
constrain the distribution of dust along the line of sight 
by applying eq. (jSJ) for various dust distributions. In the 
bottom panel of Fig. ^we show the fractional difference 
in the halo contribution to the SBD assuming half of the 
dust is uniformly distributed, and half is divided into 
two clumps at x = 0.33 and x — 0.66 (approximated as 
Gaussians with a = 0.1). The statistical errors in the 
SBD are much smaller than the fractional difference due 
to different distributions. 

4. DELAYED BURST LIGHT CURVE IN THE HALO 

In order to be able to clearly observe the delayed 
pulse of radiation from type I X-ray bursts in the X- 
ray halo, the product of the source flux and scattering 
optical depth must be large enough to provide adequate 
statistics despite smearing of the delayed burst photons 
over many kiloseconds. A burst light curve is shown by 
Thompson et al. (2005, see Fig. 1) for Ginga (GS) 1826- 
238, and can be well-approximated by an exponential 
with e-folding time <b ~ 60 s. In general, burst decays 



are shorter at higher photon energies due to a soften- 
ing of the burst spectrum during the decay, but we treat 
t\y as energy-independent for simplicity. The 60 s burst 
scale time for GS 1826-238 is relatively long; in fact, the 
majority of type I bursters have scale times less than 20 
s (see Fig. 4, Cornelisse et al. 2003). Unfortunately, the 
amount of dust towards GS 1826-238 is too low to apply 
the TS73 method, with AT H « 1.9 x 10 21 cm" 2 (Dickey & 
Lockman 1990) implying r sca ~ 0.05 at 1 keV (Predehl 
& Schmitt 1995). 

The small size of the expected signal from a type I 
X-ray burst is immediately apparent considering that 
approximately (F peak - F perBistellt )T sca i6 singly-scattered 
photons cm -2 keV -1 are expected over all halo angles 
and delay times for an X-ray burst approximated by an 
exponential. At halo angles larger than 100", the burst 
photons are smeared to such a large extent that any sig- 
nature of the delayed burst is hard to observe. To begin, 
it is straightforward to show that the time delay of pho- 
tons scattered at a fractional distance x to the source 
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We have invoked <j) ~ 0/(1 — x ), which is true for small 
angles and applies to all X-ray scattering. Although the 
time delay for a photon scattered at a known fractional 
distance and angle is simple to derive, in practice we have 
to integrate over a distribution function of time-delayed 
photons, which in the single scattering approximation is 



H(t d ,e,E) 
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t c \ + tg dQ 
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where we have defined tg = DO 2 /2c. This distribution 
can be understood by noting that the delayed light curve 
in the halo is only a function of the time delay and view- 
ing angle, although physically, the light curve is governed 
by the distribution of dust, and the differential scattering 
cross section of the dust; the latter are functions of x and 
4>, respectively. Therefore, -ff (id, 0) is simply the product 
of p{x), da((f>)/dfl, and the Jacobian determinant relat- 
ing the area elements \d(x, (f>)/d(t d , 9)\ (TS73). Although 
da/dQ is a function of the scattering angle, it is implicitly 
a function of the time delay through cf> « 0/(1 — x). An 
alternative derivation of H(td,0) is provided by Alcock 
& Hatchett (1978, see eq. [17]), who also derived formu- 
lae describing the propagation of a pulse of radiation in 
the large optical depth limit. 

The distribution of single-scattered delayed photons 
has a form similar to a decaying exponential, falling 
off more slowly at larger halo angles and for larger 
source distances. Considering that more than half of 
the space along the line of sight is covered within about 
600(L>/5 kpc)(6» h /10") 2 s (dx/dt d ~ ij 2 ), this type of 
functional behavior is not surprising. The distribution of 
delayed photons also shows that the extent of smearing 
increases dramatically with halo angle, which can also be 
seen in eq. (JSJ, where t d oc 9 2 . While the delayed light 
curve of a <5-function burst at a given angle can be ob- 
tained directly from H(td, 0), for an arbitrary light curve 
one must convolve the non-delayed light curve I(t) and 
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the distribution of time-delayed photons: 
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The delayed burst light curve (in units of counts/s) for 
singly-scattered photons for a given range of halo angles 
and energies, and for a given dust distribution and source 
distance is 

e 2 e 2 

B(t d ,AE,A9) = J J B(t d ,E,9)(S h (E)) [e r(B) - l] (11) 
e 1 e x 

xA cS (E,8)g p (E,9)2n6d8dE, 

where (Fx,b(E)) is the average burst spectrum (with the per- 
sistent spectrum subtracted out), and requiring eq. flu! to be 
normalized, which is obtained naturally if 1(f) is normalized, 
i.e. I(t) — 1 ex P (—t/th), for a type I burst approximated by 
an exponential. Figure [3] shows the normalized delayed light 
curves for X-ray bursts with ib = 60 s and E = 2 keV for 
sources at three distances with uniformly distributed dust. 
The ratio of the peak to persistent flux was chosen to be 
7, similar to that of GS 1826-238. This ratio is fairly typi- 
cal of brighter bursters (low persistent emission sources have 
much higher ratios; Cornelisse et al 2002), although the typ- 
ical range is rather wide. To facilitate interpretation of Fig. 
[3] consider the case of a unit annulus centered on 5" and a 
4 kpc source distance (top panel, dotted curve): About 12% 
of the total burst fluence would be expected in a 100 s bin 
centered at a time delay of 100 s. Small angles are shown, 
though the general trend of a smaller signal at larger angles 
due to greater smearing is clear. 

5. SIMULATIONS 

5.1. Maximizing the Signal-to-Noise Ratio 

The feasibility of using type I bursts as distance indica- 
tors was tested by simulating burst light curves in the halo, 
including both statistical and systematic errors, and finding 
the deviation of the best-fit distances. Systematic errors in- 
clude uncertainties in the calibration of the PSF, and in the 
measurements of the optical depth, the peak and persistent 
fluxes, and the scale time of the burst. Due to the small size 
of the expected delayed burst signal, it is essential to max- 
imize the S/N ratio (where the signal refers to counts due 
to the burst and the noise refers to the total Poisson noise) 
by binning the data in a favorable configuration for the sim- 
ulations. At higher energies, for example, the dust-scattered 
halo is rapidly decreasing while the PSF contribution is in- 
creasing. Not only does this binning depend on the spectrum 
and optical depth of the source, the observing telescope's ef- 
fective area and PSF, and the distribution and composition 
of interstellar dust, but also on the time after the onset of the 
burst. Due to the energy dependence of the differential scat- 
tering cross section, harder photons appear at smaller time 
delays and halo angles relative to softer photons. Our goal 
was therefore to find an effective binning of the data averaged 
throughout the dynamical time period of the delayed burst 
evolution, meaning the period when the delayed light curves 
for sources at different distances differ by more than a few per- 
cent. For uniformly distributed dust, the dynamical period is 
isolated to time delays less than ~10 ks because the delayed 
burst light curves for sources at different distances converge 
at longer time delays (see Fig. The time-averaged opti- 
mal binning was obtained by selecting an energy grid from 
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Fig. 3. — Predicted delayed burst light curves for small halo 
angles (top to bottom: 9^ = [5,10,15,20] arcsec), and for three 
source distances (dotted: 4 kpc, solid: 6 kpc, dashed: 8 kpc) as- 
suming uniformly distributed dust, 2 keV photons, and tb = 60 s. 
The j/-axis units are the percent of total burst photons appearing 
per second per unit solid angle. 



1-6 keV with 0.5 keV bins, and a radial grid from 10-160" 6 
with 10" bins. Because the cross section for X-rays scattering 
off of interstellar grains is proportional to ~_B -2 , and because 
the smearing of the burst signal increases dramatically with 
angle (t d oc 6 2 , see Fig. [3J, the lower bounds of the trial 
binning were fixed to 1 keV and 10". Using 200 s time bins, 
trial combinations of angle and energy were selected. The in- 
cluded photon energies were increased incrementally, and for 
each energy increment the included range of halo angles was 
also increased incrementally. For each trial binning, the time- 
averaged S/N ratio was calculated. We find that averaged 
throughout the dynamical portion of the delayed burst, the 
optimal annular region is about 10-50" and the optimal en- 
ergy is from about 1-5 keV for a source with r sca (l keV) ~ 1.4 
and the spectrum of GX 13+1 observed with Chandra's effec- 
tive area and PSF. This range of halo angles is not surprising 
considering the ratio plot in the lower panel of Fig. al- 
though the reader should note that we did not optimize the 
plotted ratio but rather the S/N ratio. 

5.2. Distance Measurement Method and Potential 
Accuracy 

The best-fit distance to the simulated light curves was 
measured using ^-minimization. A five-dimensional array 
composed of the predicted flux as a function of delay time 
for an arbitrary distance, halo angle, and energy was cre- 
ated. The energy and angle degrees of freedom were removed 
by integrating over the predicted single-scattered burst halo 
spectrum for GX 13+1, and the angles 10-50" assuming an 

6 Halo angles larger than 160" were not analyzed due to the 
presence of a diffuse "transfer swath" caused by the diffracted spec- 
trum, analogous to the so-called "transfer streak" centered on the 
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Fig. 4. — Top: Predicted 10-50" halo light curve and a sample 
simulated light curve after the onset of a type I X-ray burst for a 
source at 4 kpc, assuming only statistical errors. Middle: Relative 
contributions to the predicted light curve. Bottom: Histogram (0.1 
kpc bins) of the best-fit distances for one thousand simulated light 
curves, and the Gaussian fit to the distribution. 



azimuthally-symmetric halo profile. The persistent halo was 
treated as a constant. The value of % 2 was then calculated 
by comparing predicted light curves for sources at distances 
in 0.01 kpc increments to each trial of simulated data; the 
best-fit distance for each simulated light curve had the small- 
est fit statistic. One thousand trial light curves were used in 
each simulation. We first examined the potential accuracy of 
a distance measurement using a Chandra grating observation 
assuming the only uncertainties are due to Poisson noise. By 
simulating burst light curves for source distances from 3 kpc 
to 11 kpc in 0.5 kpc increments, we find that as the source 
distance becomes larger, the distribution of best-fit distances 
becomes increasingly distorted from a Gaussian due to a slight 
tail for larger source distances (forming a Gamma-like prob- 
ability density function). This makes sense because the dif- 
ferences in the delayed light curves become smaller at larger 
distances. Therefore, a negative statistical error (fewer counts 
than the average) will overestimate the distance more than 
an equal positive statistical error will underestimate the dis- 
tance. Figure [I] shows predicted and simulated 10-50" halo 
light curves, the relative contributions of the halo and the 
PSF, and a histogram of the best-fit distances for a source 
distance of 4 kpc. We find the fractional error on the in- 
ferred distance to a source 4 kpc away is j^% for a single 
burst (68% confidence). At 6 kpc and 8 kpc, the errors are 



% and _ 21 % respectively. The positive and negative errors 



+29 c 
-19' 

increase approximately linearly with distance. For a source 
with a smaller optical depth to scattering than GX 13+1, e.g. 
T sca (l keV) ~ 0.5, roughly five aligned and summed bursts 



would be required to obtain similar accuracy. By including 
realistic systematic uncertainties of 10% for the scattering 
optical depth, 20% for the PSF contribution, and 5% on the 
flux and scale time of the burst, the potential accuracy de- 
creases by only ~10% for distances between 3 kpc and 11 
kpc. Future large area observatories with adequate spatial 
resolution (e.g., XEUS and Constellation-X) will effectively 
minimize the Poisson noise, in which case the distance mea- 
surement accuracy is limited by systematics. By simulating 
data with systematic errors only (which is also equivalent to 
stacking many bursts), the uncertainty in the distance mea- 
surement is 7% at 3 kpc, increasing approximately linearly 
with distance to 12% at 11 kpc . 

The simulations above assumed a constant persistent halo, 
but due to the nature of the accretion processes in LMXBs, 
there is likely to be source variability on all time scales. With 
reliable flux measurements (from a grating observation), han- 
dling the PSF is straightforward because its contribution is 
essentially a constant fraction of the point source flux (see eq. 
|2|), and there is no time delay associated with it. Likewise, 
halo variability can be handled by including additional terms 
in the light curve I(t) in eq. IjlOI . and appropriately chang- 
ing the normalization. As a worst case scenario, however, 
assume that variability exists on 200 s scales, yet we still as- 
sume the persistent halo is constant in the calculation of the 
delayed light curve. To test the effect of this assumption, 
we calculated delayed light curves for halo angles of 10" to 
50" in 10" increments for time-averaged constant source light 
curves (on a scale much longer than 200 s) but with stochastic 
variability on 200 s scales with rms amplitudes of 10%, 20%, 
30%, and 40%. Because this variability is spread out over 
many kiloseconds due to scattering along the line of sight, 
the rms amplitude is smaller in the halo. Interestingly, the 
halo rms amplitudes seem to be only weakly dependent on the 
halo angle, with corresponding rms amplitudes are ~2%, 4%, 
5-7%, and 6-9% respectively. These percentages were added 
to the statistical and systematic errors for the persistent halo 
contribution in our simulations. For the worst case of 40% 
point source variability, we find the inferred distance errors 
increases by about 15% for a grating observation, e.g., for a 6 
kpc source distance and a simulation including systematics, 
the error increases from about 30% to 35%. This can also 
account for dust distributions that are overall uniform, but 
patchy on scales i < 1. 

A non-grating observation with Chandra increases the ef- 
fective area and signal by about a factor of 4 (see Fig. 2), but 
the reduction in Poisson noise is offset by larger systematic 
uncertainties in the flux measurements and the burst scale 
time. Without the dispersed grating counts, one is required 
to extract spectral and flux characteristics from the source 
transfer streak (e.g., Smith et al. 2002, Clark 2004), gen- 
erally a much coarser and poorly calibrated method. More 
importantly perhaps, a non-grating observation would signif- 
icantly increase the extent of pile-up, forcing one to look at 
larger halo angles where the evolution of the delayed burst 
occurs on longer time scales. 

6. CONCLUSIONS 

In this paper, we studied the applicability of the TS73 
method to the delayed signature of type I X-ray bursts in the 
halo. The signal produced by an X-ray burst in the halo is in- 
variably small. Regular type I bursts last ~5-150 s, but the 
scattering geometry, which governs the distribution of time 
delays, assures that these photons will be spread out over 
time delays greater than ~10 ks for uniformly distributed 
dust. Nevertheless, if the physical parameters governing the 
light curve of the burst in the halo are well-known, the statis- 
tics provided by Chandra's effective area can yield a distance 
measurement to about 25% accuracy with a single burst. By 
stacking many bursts, systematic errors likely limit the poten- 
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tial accuracy to roughly 10%. Although non-grating obser- 
vations have improved statistics due to larger effective area, 
pile-up precludes accurate flux measurements (although the 
transfer streak can be used) and require the use of larger halo 
angles, decreasing the S/N as the evolution of the delayed 
burst light curve occurs on longer time scales. Source vari- 
ability on 200 s scales with 40% rms amplitude increase the 
distance uncertainty by ~15%. The quoted errors assume 
very favorable source characteristics: (1) a bright source with 
a flux ~0.6 photons/cm 2 /s during persistent emission from 
1-5 keV, (2) a large optical depth to scattering (r sca ~1.4 at 
1 keV), both of which increase the S/N ratio of the delayed 
burst light curve in the halo, (3) a burst with large total flu- 
ence, characterized by t b = 60 s and FpeakZ-Fpcrsistcnt = 7. In 
the general case, the uncertainty in the inferred distance will 
be somewhat larger. Finally, note that the interstellar grain 
model represents an additional source of systematic error. 

In many cases the distance to Galactic X-ray sources has 
been estimated using other methods, e.g., using the require- 
ment that the inferred peak luminosity is sub-Eddington, or 
using optical photometry of the (low-mass) binary compan- 
ion, and the potential accuracy of the TS73 method applied 
to bursters may only provide additional distance constraints 
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